The Lifshitz Transition in d-wave Superconductors 
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The BCS to BEC evolution has been recently the focus of studies in superconductors and cold 
atomic gases. For a d-wave system, we show that a Lifshitz transition occurs at a critical particle 
density which separates two topologically distinct phases according to their quasiparticle excitation 
energies: a BCS-like gapless superconductor in the higher density limit and a BEC-like fully gapped 
superconductor in the lower density limit. This transition is second order according to Ehrenfest's 
classification, but it occurs without a change in the symmetry of the order parameter, and thus can 
not be classified under Landau's scheme. To illustrate the nature of the transition, we compute the 
compressibility and the superfluid density as functions of particle density. 
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I. INTRODUCTION 

The evolution from BCS to Bose-Einstein conden- 
sation (BEC) superconductivity/superfluidity has at- 
tracted considerable amount of interest in the condensed 
matteri^i^iMi^ and atomic physics^'^'® communities. In 
the atomic physics community, this interest resulted from 
the possibihty of studying condensation phenomena in 
fermionic atomic gases, where the scattering length (and, 
thus, the effective interaction strength) can be tunned 
via Feshbach resonances for a given densitySi^. In the 
condensed matter physics community, the interest in the 
BCS to BEC evolution arouse in the context of high-Tc 
superconductivity'^, where the nature of the supercon- 
ducting and normal states changes as function of car- 
rier concentration. Furthermore, experimentalists seem 
to start having some control over the carrier concentra- 
tion using electronic doping via ferroelectric oxidesSi. The 
fundamental issue that needs to be addressed is whether 
there is a quantum phase transition in the evolution from 
a BCS to BEC ground state, as particle density or scat- 
tering length (interaction strength) are varied. In this pa- 
per we show that for a two-dimensional d-wave supercon- 
ductor there is a Lifshitz transition between the BCS and 
BEC ground states. The transition is second order ac- 
cording to Ehrenfest's classification of phase transitions, 
but it occurs without change in symmetry as Landau's 
classification would demand. The quantum phase transi- 
tion found has a logarithmic singularity in the compress- 
ibility at a critical concentration for fixed interaction. 
Furthermore, near the critical point the low temperature 
superfluid density exhibits a dramatic change in behav- 
ior. 

This manuscript is organized as follows. In Section 
II, we discuss the Hamiltonian and the interparticle po- 
tential used in our model. In Section III, our functional 
integral calculation is presented, and the order parame- 
ter and number equations are derived at the mean-field 
level. The analysis of Gaussian fluctuations about the 
saddle point solution is also performed in this section. In 
Section IV, our results for the electronic compressibility 
and the phase diagram of the system are shown, while 



in Section V an analogy with the Lifshitz transition in 
metals is discussed. The superfluid density is analyzed 
in Section VI and, finally, our concluding remarks are 
summarized in Section VII. 



II. HAMILTONIAN AND INTERACTION 
POTENTIAL 



We study a two-dimensional continuum model of 
fermions of mass m described by the Hamiltonian (h = 

kB = 1) 



(1) 



k,(T 



k,k',q 



where &kq = V'-k+q/2T'0k+q/2i and ek = k"^ /2m. We 
consider the following separable potential in k-space. 



kk' 



-Adr(k)r(k'), 



(2) 



which includes only the dominant angular momentum 
channel, assumed to be d-wave. The interaction term 
can be written as r(k) ~ h{k)g{'k), where h{k) — 
(fc//ci)^/[l -|- k/ko]^/^ controls the range of the interac- 
tion, 5(k) — cos(2(^) sets its angular dependence, and 
ip is the momentum angle in polar coordinates^". In 
this case, feg ^ Rq^^ where Rq plays the role of the 
interaction range, and both fco and ki set the momen- 
tum scales in the short and long wavelentgh limits. We 
work under the assumption that the system is dilute 
enough, i.e., kp^^^ <C fcg. When computing physical 
properties throughout the manuscript, we scale momenta 
by fci?^^^, energies by ei?,^^^ = kp^^^/2m, velocities by 
"'^i^max — ^^'J^max/'^: ^ud partlclc dcuslty n by nmax/27r, 
— kp^^^/2TT. If we choose, for instance, 

{Ro « 0.32A) and define kp^^^ = fco/10, 



where n^ax 

ko = VioA- 

then rima^ r 



1.59 X lOi-^cm 



We next discuss the ef- 
fective action and analyze the effects of Gaussian fluctu- 
ations about the saddle point solution. 
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III. EFFECTIVE ACTION AND GAUSSIAN 
FLUCTUATIONS 

The partition function Z at a temperature T = 
f3~^ is written as an imaginary-time functional integral 

with action S = /q^ (^TEk,<T V'L('^)(^^ ~ ^J')i'^c{T) + 
Ti]. Introducing the usual Hubbard-Stratonovich field 
(/>q(r), which couples to V^V'^i E^nd integrating out 
the fermionic degrees of freedom, we obtain Z = 
/ exp(— S'cff [0, with the effective action 

given by 



'S'cff = 



dT 



U{r) 



E 

k.k' 



k,k 



where = Ck - U{t) = J^i, l0k(r)|V^> and Gkj^,(T) 
is the (inverse) Nambu propagator, 



Gkl'(^) 



-(9r +Ck)'5k,k' Ak,k'(T) 



(3) 



with Ak,k'(T) =0k-k'(r)r((k + k')/2). 

To study the BCS to BEC evolution it is necessary 
to go beyond the standard saddle point approximation, 
and include at least Gaussian fluctuationa^. Assuming 
(/)q(T) = Apf^q^o + '7q(''")j ^nd performing an expansion in 
^eff to quadratic order in 77, one obtains 



^Gauss = So[Ao] + i ^^t(^)M(g)7y(<j), 



(4) 



where So is the saddle point action, 77^ = [ry*(g), r7(— g)], 
and q = {q,iqm) with iq^ — i2rm:/(3. The inverse fluc- 
tuation propagator M is a 2 x 2 matrix to be defined 
later. 

The Saddle Point Equation: Starting with the sta- 



tionarity condition [6Sos/S4'qi'''')]i 



0, Fourier trans- 



forming from imaginary time to Matsubara frequency, 
ikn — i{2n l)7r//3, and performing the frequency sum, 
the saddle point condition can be written as 



1 



V^^Mtanh^'^^'^ 

4" 2£;k 



(5) 



where i?k = -y/^^ -I- J\'^T'^{'k) is the quasiparticle excita- 
tion energy, with ^k = ek ^ Fig. ^ shows a plot of 
±£^(k, /i) as a function of k for three different values of 
the chemical potential ^i. Notice that the Dirac cones 
collapse at = 0. 

The Number Equation: Using the thermodynamic re- 
lation N = —dfl/dfi, and taking f^causs = fio + ^^fluct, 
where flo = S'o[Ao]//3 and f^fluct = /J-^ lndet[M(q)], 
one can write the number equation as 



A^Gauss = A^O + A^fluct, 



(6) 



where A^o — —drio/dfi — 2^j^nk, with rik ~ 
1/2 [1 — ^ktanh(/3i?k/2)/£'k] being the momentum dis- 
tribution, and 



-9(detM)/9/i 



detM(q,ig„) 



(7) 



being the fluctuation contribution to Nqi^^ss- At T — 0, 
we find well defined Goldstone modes for all couplings, 
provided that |q| is sufficiently small. This collective 
mode appears as a pole in the two-particle excitation 
spectrum determined by det[M(k, z)] = 0, with ig„ — > z. 
This pole has the form z = c|q| — zd|qp, where c > is 
the speed of sound and d > is the damping coefficient. 
For fi < 0, d vanishes and the contribution from the 
collective mode pole z — c|q| dominates at sufficiently 
low temperatures, since the two-particle excitations are 
gapped. For yu > 0, cZ becomes positive, and the spectrum 
of two-particle excitations is gapless due to the presence 
of the Dirac points. Thus, unlike the s-wave case. Lan- 
dau damping occurs even at T = 0. However, for small 
|q|, the Goldstone mode is underdamped, i.e., the real 
part of the pole dominates. Therefore, for either /i > 
or /It < 0, these fluctuation effects lead to 



A^fiu 



C(3) 1 dc 



rp6 



N 2 ' 

max ^ 

which vanishes in the limit of T — )■ 0. Thus, we recover 
the rf-wave equivalent of Lcggett's variational results^-^ 
at T = for the saddle point and number equations in 
the context of '^He. 

Denoting ^k and by ^ and i?, and ^k+q and -Bk+q 
by ^' and E' , the matrix elements of M become 



2EE' 



E+^E' 



(£+£')r^(k+q/2) 

(iq„)^-iE+E')^ 



2EE' 



r"(k+q/2) 

(ig,.)^-(-E+B')^ 



iqn 



M?^{q)^-Eu{^ 

M (n\--\^ / Agr^(k+q/2) \ / (i;+g')r(k)r(k+q) 
iviuyq) — l^y. y J y (iq„Y-(E+E'Y 

where the element Mii(g) = M^^{q) + Mg(g) was split 
into an even [E) and an odd (O) part in iq^ and 
^^21 (9) = Mi2{q) and M22{q) = Mii(-q). Performing 
the analytic continuation iqn ^ a; -I- lO"*" and defining 
rj{q) — \rj{q)\e.yi\){i(j){q)) — {\{q) + i9{q))/^/2, the correc- 
tion to the effective action assumes the form 



(A* 6*) 



+ M12 



Affi - M12 




(9) 



where \{q) and 6{q) are real functions that may be identi- 
fied with amplitude and phase fluctuations, respectively. 
Since the excitation spectrum is gapped for /i < 0, one 
can make a direct small q and uj expansion, result- 
ing in + Afi2 ^ A - Bu^ + Q|q|2 + . . . , = 
Cuj + . . . , a.nd M^^ - M12 = -Du'^ + R\q\'^ + ..., where 
the coefficients are given by A = X^k T^Ck) A^/2£'3, 

B = j:^2er^{k)/i6E^ c = Eke^2(k)/4i^;^ d = 

Ek2r2(k)/16i;3, etc. 



(al (b) ' (c) 

FIG. 1: Plot of ±ii'(k, /i) for (a) > 0, (b) ^ = and (c) /i < 0. Notice the collapse of the Dirac cones at = 0. 



IV. ELECTRONIC COMPRESSIBILITY AND 
PHASE DIAGRAM 

We now discuss the behavior of the electronic com- 
pressibility and phase diagram of the system as the /i = 
point is crossed. A simple analysis of i?k indicates that, 
for /i > 0, i?k is gapless, while for fi < 0, is gapped. 
Simultaneously, there is a massive rearrangement of the 
momentum distribution rik as /i passes through fic — 0, 
leading to the vanishing of the first derivative of /i with 
respect to the density n = N/L"^ at n — ric- An impor- 
tant thermodynamic quantity that depends directly on 
both rik and i?k is the isothermal electronic compressibil- 
ity K, defined by n'^K = [—d^Q/dfi"^] /L"^. This quantity 
can be written as 

n^K = ao + afluct, (10) 

where ao — dno/dfj, and aauct — dufiuct/dfi. The 
term ao at T = can be explicitly rewritten as ao — 
(4/i^)Ek"k(l - nu.)/Ek, while the term afluct oc 
vanishes as T — > 0. In the vicinity of /x = /Xc = 0, k 
diverges logarithmically at T = as 

K « [-ci In |1 - ra/nd + C2] (11) 

in the d-wave case, suggesting the existence of a quan- 
tum phase transition (QPT). This singular behavior of 
the compressibility is shown in Fig. 12 and the corre- 
sponding phase diagram is shown in Fig. 13 together with 
a plot of /i as function of n (inset). The critical line 
fi — fic ^ in the Ad X 71 space separates a BCS-like 
region, where Dirac points exist and there are gapless 
excitations {p, > fic), from a BEC-like region which is 
fully gapped (/x < /Zc). This critical line ends at the 
two-body bound state threshold X*^ = S (see Fig. 3) 
when n ^ 0. Thus, the QPT between the BEC and 
BCS regimes requires that \d > AjJ. In the s-wave case 
there is no QPT, since k is continuous and changes in 
nk are always smooth. Notice that the contribution of 
the collective modes to n and to k at T = vanishes 
identically. Thus, Gaussian fluctuation effects of the su- 
perconducting order parameter are not important for the 
T = electronic compressibility. This divergence in the 
second derivative of O at T = signals a second order 
quantum phase transition, according to Ehrenfest's clas- 
sification. However, the symmetry of the order parameter 




0.5 1 1.5 



n 

FIG. 2: Plot of n^K (in units of nmax/27reF,nax) ^-nd the elec- 
tronic compressibility k (in units of 27r/nmaxeF,„ax) ^ func- 
tions of particle density (in units of nmax/27r) for fixed inter- 
action strength and ko — ki = 10fcF„,ax- 

does not change, and a Landau symmetry classification of 
the phase transition is not possible. So, if the symmetry 
is not changing at the transition point, what is? 

V. THE LIFSHITZ TRANSITION 

To answer this question we make an immediate connec- 
tion to the Lifshitz transition— in the context of ordinary 
metals at T = and high pressure. The Lifthitz tran- 
sition should not be confused with the Lifshitz point^'i', 
where a finite temperature phase transition occurs sep- 
arating the high temperature disordered phase, the spa- 
tially uniform ordered phase, and the spatially modu- 
lated ordered phase. In the conventional Lifhitz transi- 
tion, the Fermi surface e(k, P) = Ep changes its topology 
as the pressure P is changed. For an isotropic pressure 
P, the deviation AP — P — Pc from the critical pressure 
Pc is proportional to A// = /j, — jic- A typical example 
of the Lifshitz transition is the disruption of a neck of 
the Fermi surface which leads to a non-analytic behavior 
of the number of states N{fi) inside the Fermi surface. 
In this case, N{^) behaves as A(/Ltc) + — fXcl^/"^ for 
/i < /Zc, and as A{pc) for fj, > ^c, in the vicinity of ^c- 
Here, n = (3/2)P|/i - /id^/Vn^, where Uc = Nc/V. No- 
tice that K is non-analytic, but it is not singular. The 
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FIG. 3: Phase diagram of particle density n (in units of 
/27r) versus interaction strength Ad (in units of g^^^, 
where g2u is the density of states in two dimensions) for 
ko = ki — 10feF„ax- The solid line {/i = 0) separates a gapless 
regime (/i > 0) from a fully gapped regime (^ < 0). Inset: 
Chemical potential jj, (in units of ep^^^) as function of particle 
density n across the transition region. Notice that jj, vanishes 
at n = Tic = 1. 



quantity that signals a phase transition in this case is 
not K, but the thermopower Q, which is proportional 
to —d\n{n)/dfj,, thus leading to Q c>c — |A/i|~^/^. In the 
conventional Lifshitz transition, the system lowers its en- 
ergy by AE cx — |A^|^/^ cx — |AP|^/^, and the transition 
is said to be of second-and-half order—. 

For the Lifshitz transition in d-wave superconductors, 
there is already a non-analytic and singular behavior in 
dn/dfi, and thus in k. The singularity in n is logarith- 
mic as indicated above, since the system lowers its en- 
ergy by AE oc — I /i — /icp In I /i — /ic I- This logarithmic 
contribution originates from the simultaneous collapse of 
the four Dirac points at k = 0, which produces a gap in 
the excitation spectrum i?k and a massive discontinuous 
rearrangement of the momentum distribution rik m the 
ground state as /i — > = 0. A direct topological anal- 
ogy with the standard Lifshitz transition can be made by 
noticing the collapse of the Dirac cones ai fj, = fic (and 
the disruption of a neck for n < ^c) m the excitation 
spectrum of the system, as shown in Fig. ^ • 



VI. SUPERFLUID DENSITY 

Now we analyze the behavior of the superfluid density 
tensor pij (T, n) across the Lifshitz transition, given by 



(12) 



where rik is the momentum distribution and Ik = 
(2T)-isech^ (£'k/2r) is the Yoshida distribution. No- 
tice that pxx = Pyy = P, while p^y = Pyx = 0. 
Due to Galilean invariance of our continuum model, 
p{T = 0) n/m is well-behaved as the critical point 
is crossed. However, dp{T ~ 0)/dp = n^K/m diverges 



■a 




FIG. 4: Plot of the superfluid density zero temperature slope 
(in units of nmax/27reF„ax) a function of n (in units of 
?imax/27r). Inset: Ap (in units of nmax/27r) as a function of T 
(in units of eF„,ax) for various particle densities. 



at n — ric- Using our energy and momentum scales, we 
define Ap = {mp{T) — n)/2-K. In Fig. ^ we show Ap for 
various concentrations (inset), and the zero temperature 
slope {m / 2T:)[dp / dT]T=Q as a function of particle density. 
Notice the discontinuity in the zero temperature slope as 
the critical point n = Uc 'is approached. 

We have also calculated Ap analytically at low tem- 
peratures in the case of very short range interactions 
{ko oo). In the BCS limit, Ap = -2Tfc2/Ao(n), due 
to the nodal structure of the c?-wave symmetry. At the 
critical point ^ = (n = nc), Ap = — \Ti{2)F{rj)T , where 
F{ri) = {1 + ^?^)"^/^ with 7? = Ao{n)/kl. In the BEG 
limit, Ap = -(8/7r)exp(-|//|/r)F(?7). This exponential 
behavior reflects the appearance of a full gap in the ex- 
citation spectrum for n < Uc- 

For all values of p, there is further reduction of p{T) at 
low T due to Goldstone modes leading to Apc = —AT^, 
which obviously does not affect the T = slope of Ap 
through n — Tic- 



VII. SUMMARY 

We have proposed the existence of a Lifshitz transi- 
tion at T = in clean c?-wave superconductors, where a 
topological change in momentum space is responsible for 
a non-analytic behavior of the electronic compressibil- 
ity and of the zero temperature slope of the superfluid 
density. We conclude by suggesting that the search for 
this transition may now be possible using the so-called 
ferroelectric field effect transistor (FFET)Si, where some 
control over the particle density in cuprate superconduc- 
tors may be achieved without chemical doping. In ad- 
dition, it may be possible to investigate the occurrence 
of this transition by measuring the superfluid density as 
a function of dopingli, or through a direct measurement 
of the electronic compressibility as a function of parti- 
cle density, as was done in the study of metal-insulator 
transitions—. 
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